Cosmology from MAXIMA-1, BOOMERANG & COBE/DMR CMB Observations. 
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Recent results from BOOMERANG-98 and MAXIMA-1, taken together with COBE-DMR, 
provide consistent and high signal-to-noise measurements of the CMB power spectrum at spher- 
ical harmonic multipole bands over 2 < ^ < 800. Analysis of the combined data yields 68% 
(95%) confidence limits on the total density, f2tot — 1-11 ± 0.07 (to }2)' ^^e baryon density, 
^bh"^ ~ 0.032t^;E;|^^ (taoos). and the scalar spectral tih, n, ~ l.Olto.o? (tall)- These data 
are consistent with inflationary initial conditions for structure formation. Taken together with 
other cosmological observations, they imply the existence of both non-baryonic dark matter and 
dark energy in the universe. 

PACS numbers: 98.80.-k; 98.70.Vc; 98.80.Es; 95.85.Bh 



Measurements of the angular power spectrum, Cg, of 
the Cosmic Microwave Background (CMB) have long 
been expected to enable precise determinations of cos- 
mological parameters Q . The CMB power spectrum de- 
pends on these parameters, as well as the scenario for 
the generation and growth of density fluctuations in the 
early universe. Evidence for structure in the CMB of the 
character predicted by adiabatic inflationary models has 
been mounting for the past decade and was convincingly 
detected in 1999 [|, |, |. The recent BOOMERANG-98 
(B98) D and MAXIMA-1 [§ CMB anisotropy data pro- 



vide a significant improvement in the determination of Ci. 
This letter jointly analyzes these two data sets, incorpo- 
rating COBE-DMR 1^] and other cosmological informa- 
tion to obtain further estimates of several cosmological 
parameters. 

MAXIMA-1 and B98 have produced independent 
power spectra from patches of sky roughly 90° apart, on 
opposite sides of the galactic plane (||, ||] and references 
therein) . These data provide the first narrow-band detec- 
tions of the power spectrum from 400 ^ £ ^ 800, where 
further acoustic peaks are expected in adiabatic models. 
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Each spectrum shows a well-defined peak at multipole 
I ~ 200, followed by a relatively flat region extending to 
the highest multipoles reported (Figure]^, top). These re- 
sults have been interpreted as supporting the inflationary 
theory of structure formation with adiabatic initial con- 
ditions, and allow the first precise CMB measurements 
of other parameters such as the baryon density ^, ^ . 

Comparison & Calibration. The B98 and MAXIMA- 
1 angular power spectra |l2j are shown in the top panel 
of Fig. |l|, along with some best-fit models. The B98 
(MAXIMA-1) data cover I = 25-625 (36-785) with a 
resolution 5t — 50 (75). The DMR data provide informa- 
tion at low Each data set has a calibration uncertainty 
(20% for B98, 8% for MAXIMA-1, Icr in C() and a beam 
uncertainty (10% for B98, 5% for MAXIMA-1) that are 
not included in the errors plotted in Fig. |l|. We take 
these uncertainties into account by allowing the overall 
amplitude (calibration) as well as an ^-dependent am- 
plitude (beam) to vary for each spectrum, weighted by 
the distribution of errors around the nominal beam and 
calibration values. 

We have combined the two data sets in the bottom 
panel of Fig. ^ We approximate the likelihood in the 
manner of ll| , treating the individual bandpowers as 
the parameters to be determined. The combined power 
spectrum requires an increase in the calibration by a fac- 
tor 1.14 ± 0.10 for B98 and a decrease by 0.98 ± 0.08 for 
MAXIMA-1, along with a beam rescahng of 1.07 ± 0.09 
for B98 and 0.99 ± 0.05 for MAXIMA-1. 

We define a goodness-of-fit parameter, = — 21n£, 
in terms of the likelihood function, C, which reduces to 
the usual for a Gaussian. The data shown in the top 
panel of Fig. |l| have x^/DOF = 11.6/8 with respect to the 
combined spectrum of the bottom panel (22 data points, 
10 bandpowers and 2 beam rescalings and 2 calibration 
parameters give 8 degrees of freedom [DOF]) when DMR 
is not included, and a x^/DOF = 34/29 with DMR. This 
indicates that the results are each consistent with a single 
underlying power spectrum. 

Cosmological Parameters. We consider a subset of 
parameters describing a Friedmann-Robertson- Walker 
Universe with adiabatic initial conditions. The present- 
day density pi of a component i is fli = SirGpi/ {3Hq) 
where Hq — lOO/i km/s/Mpc is the Hubble constant; 
thus Uih^ is a physical density, independent of Hq. We 
consider the density of baryons, fit; the density of cold 
dark matter (CDM), ilc, the total density of matter, 
ri„i = ilc + flb; the effective density of a cosmological con- 
stant, Q\; and finally the total density iltot = + ^A, 
ignoring other possible components. The initial spec- 
trum of density perturbations is described by an ampli- 
tude factor, Cio, multiplying the CMB spectrum, and the 
spectral tilt of scalar (density) perturbations, rig (defined 
so the initial three-dimensional perturbation spectrum is 
P{k) oc fc"=). We also consider the optical depth to the 
epoch of reionization, tc. Our parameter space is thus 



MAXIMA-1 o 

Boomerang — 98 x 
DMR ♦ 




400 

multipole / 

FIG. 1; CMB power spectra, Ce = £{£ + l)Ce/2ir. Top: 
MAXIMA-1, B98 and COBE-DMR. Bottom: maximum- 
likelihood fit to the power in bands for the three spectra, 
marginalized over beam and calibration uncertainty. In both 
panels the curves show the best fit model in the joint parame- 
ter estimation with weak priors and the best fit with fitot ~ 1. 
These models have {fitot, f^A, ftbh^, ilch^ ,ns,Tc}~{1.2, 0.5, 
0.03, 0.12, 0.95, 0}, {1, 0.7, 0.03, 0.17, 0.975, 0}. They re- 
main the best fits when the large scale structure prior |2^] 
is added, and when the SN prior [^2| is added the f2tot = 1 
model becomes the best fit in both cases. 



{fltot, ^bh^ n,h^, n„ TC, InCio} 

Because many of these parameters affect the spectrum 
in highly correlated (in some cases almost degenerate) 
ways, limits on any one parameter are necessarily a func- 
tion of the constraints, implicit and explicit, that one as- 
sumes for the other parameters. Thus, all such limits 
must be understood in the context of the specific ques- 
tion that one asks of the data. 

Inflation predicts fitot = 1- Fig. |l| shows that the data 
are consistent with this prediction: there exist fltot = 1 
models that fit the data well with quite reasonable choices 
for all of the other parameters. Simultaneously max- 
imizing the likelihood for all parameters at once gives 
i^tot = 1.2. As is clear from the Figure, the global best- 
fit model and the best-fit flat model both fit the data with 
X^/DOF < 1. Extracting quantitative conclusions from a 
simple x^ analysis is complicated by the fact that param- 
eters are correlated with one another and are bounded 
by other cosmological observations. This issue will be 
explored elsewhere. 

In this paper, we constrain parameters individually 
by marginalizing (integrating the posterior distribution) 
over all other parameters (including those describing un- 
certainty in the beam and calibration). We apply the 
techniques used on B97 [|, B98 § and MAXIMA-1 |] 
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to the combined data. We calculate the likelihood of a 
suite of models given the DMR, MAXIMA- 1, and B98 
results, taking into account the correlations within each 
data set, as well as the non-Gaussianity of the likelihood 
functions 0, We use the median, 16% and 84% in- 
tegrals as central and ilcr error estimates respectively. 
95% limits are approximately twice the quoted errors. 

The most likely value of each parameter calculated in 
this way will, in general, be different from that found 
by multi-dimensional maximization. Marginalization has 
the advantage of taking into account the likelihood struc- 
ture over the entire space, weighted by the likelihood 
value and the volume of the parameter space. It has the 
disadvantage of allowing models whose validity we may 
doubt to influence the final result. Hence, it is crucial 
to restrict the range that each parameter may take on 
through the prudent use of "prior" constraints, and to 
test how the results change under a change of priors. 

The finite range of each parameter that is included in 
our model database acts as a uniform (tophat) prior prob- 
ability density [|l^. We further restrict the analysis to 
the cosmologically interesting regime of 0.45 < h < 0.90, 
age to > 10 Gyr, and > 0.1; we refer to this com- 
bination as our "weak prior." Without such restrictions 
the data sets allow pathological low-sound-speed mod- 
els with strong positive curvature {flk < 0), very high 
baryon density, and very young ages. Parameter degen- 
eracies allow such models to contribute to or dominate 
the likelihood j|, |l^. The results presented here were 
checked using the methods of |^ and found to be in 
good agreement ]l8| . 

Results are shown in Table ||. The combined data 
yield parameters consistent with those derived individ- 
ually: the curvature is close to flat, the spectral in- 
dex is close to unity, and the baryon density is about 
2(7 above (with respect to our CMB confidence limits) 
the favored Big-Bang nucleosynthesis (BBN) value of 
ribh'^ = 0.019 ± 0.002 This high value suppresses the 
second acoustic peak relative to standard CDM models. 

In Figure || we plot likelihoods for several parameters 
using the two data sets independently and for the 
joint analysis. We also plot likelihood contours, showing 
how the CMB primarily restricts fitot = + ^A- Pro- 
jecting this figure onto either axis gives the likelihood for 
the quantities separately. This indicates the danger of 
individual parameter estimates in the presence of strong 
correlations: for example, neither Qm nor Q\ are well- 
determined individually by the CMB data alone. 

We have also investigated the effects of applying vari- 
ous prior probabilities Overall, the results from the 
combined data are somewhat less dependent on the pri- 
ors than cither individually. We find that any prior re- 
striction that docs not seriously contradict the locus of 
good-fitting models has very little effect on the parame- 
ter estimates. Incorporating previous CMB data ^, |l^ 
has very little effect. Restricting the Hubble constant to 



TABLE I: Parameter estimates from the two data sets [ |l8[ , 
and tlie combined data, using the weak prior (0.45 < h < 
0.90, to > 10 Gyr, flm > 0.1). Below the line, we restrict 
the parameter space to fitot = 1 and add other cosmological 
information. Central values and la limits are found from the 
50%, 16% and 84% integrals of the marginalized likelihood. 
Tc and Q,A are not constrained by the data. 
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0.71 ± 0.08 |2C(] does push the curvature closer to 
fltot = 1-O5^o4- A strong BBN prior (fife/i^ = 
0.019 ± 0.002 |l9|j, disfavored by the current CMB 
data, pushes other quantities to compensate, such as 
n, = 0.89 ±0.06. 

We determine parameters which are functions of those 
used to explicitly define our parameter space by calcu- 
lating their means and variances over the full proba- 
bility distribution. We find h = 0.57 ±0.11 and age 
to — 14.6 ± 2.0 Gyr. Restricting to fitot = 1 gives 
h = 0.75 ± 0.10, to = 11.8 ± 0.8 Gyr. To illustrate 
the effect of the priors alone, we give the following 
limits without including CMB data: h = 0.63 ± 0.11, 
to = 12.3 ± 1.9 Gyr. 

The combined data improve the constraint on the mat- 
ter density and cosmological constant, but we caution 
that these results are influenced by the weak prior alone 
(e.g., Qch'^ = 0.18io;io without CMB data). However, 
prior information from other cosmological data tight- 
ens these limits considerably. We consider the con- 
straints on the power spectrum from observations of 
large-scale-structure (LSS) |2^ and from observations 
of distant Supernova la (SNIa) |Q as in the bot- 
tom of Fig. |. CMB+LSS gives (r2™,r2A) = (0.49 ± 



0.13,0.63 



+0.0K;. 

-0.097) 



(note that flm is a function of the ex- 



plicit database parameters, and thus its confldence lim- 
its are calculated as discussed above) CMB+SNIa gives 
(0.35 ± 0.07, 0.75lE]-g?); combining all three gives (0.37 ± 
0.07, 0.71 ± 0.05 for' CMB+LSS+SNIa. We note that the 
combination of CMB±LSS is about as restrictive as — and 
compatible with — the combination of CMB+SNIa. 

Conclusions. The MAXIMA- 1 and B98 data are con- 
sistent with one another over the range of overlapping 
coverage in £. The consistency between the two data sets, 
obtained by different experiments using different obser- 
vation strategies on different parts of the sky, eliminates 
many sources of systematic error as cause for concern. 

These data, together with those of COBE-DMR, sup- 
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FIG. 2: Likelihood functions calculated using the weak 
prior. Tbp: likelihoods from DMR+B98, DMR+MAXIMA- 
l(M-l) DMR+MAXIMA-1+B98 and DMR+MAXIMA- 
1-I-B98+LSS (LSS is the large-scale structure prior (il)). Bot- 
tom: the likelihood in (Sim, ^^a)- Shaded contours nearly par- 
allel to flm+^A ~ 1 are one-, two-, and three-sigma limits, de- 
fined as the equivalent likelihood ratio for a two-dimensional 
Gaussian distribution, from DMR-I-B98-I-MAXIMA-1 with 
weak priors (left) and DMR-hB98-hMAXIMA-l-hLSS (right). 
Contours labeled "SNIa" are from high-redshift supernovae 
observations , and the final heavy set of contours are con- 
straints from the product of the two distributions. 



port the chief predictions of the inflation paradigm, that 
the geometry of the universe is flat, and that the ini- 
tial idensity perturbations are scale- invariant, and with 



and the theory of BBN ||. 

The data analyzed here present only a suggestion of the 
expected second acoustic peak in the CMB power spec- 
trum. A detailed mapping of that £ region with higher 
signal-to-noise may yield more convincing detections of 
acoustic oscillations, providing yet more evidence for the 
adiabatic inflationary paradigm and measurements of the 
cosmological parameters. Already, the combined data 
begin to limit the CDM density, driven by the B98 con- 
straints on the first peak and the MAXIMA- 1 constraints 
over 650 ^ ^ ^ 800. Both teams are analyzing addi- 
tional data which may significantly reduce the errors in 
the region of the power spectrum where further peaks are 
expected. 
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